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Abstract 

Higher-dimensional braneworld models which contain both bulk and brane 
curvature terms in the action admit cosmological singularities of rather un- 
usual form and nature. These 'quiescent' singularities, which can occur both 
during the contracting as well as the expanding phase, are characterised by 
the fact that while the matter density and Hubble parameter remain finite, 
all higher derivatives of the scale factor (a, a etc.) diverge as the cosmological 
singularity is approached. The singularities are the result of the embedding 
of the (3+l)-dimensional brane in the bulk and can exist even in an empty 
homogeneous and isotropic (FRW) universe. The possibility that the present 
universe may expand into a singular state is discussed. 



In this letter, we would like to draw attention to the fact that braneworld theory admits 
cosmological singularities of very unusual form and nature. The theory we work with in this 
paper is described by the action 

S = m 3 J2\[ (ft-2A0-2/ k} + [ (m 2 R-2a)+[ L(h a/3 ,<f)) , (1) 

■ L^bulk J brane J J brane ^ ' J brane 

the notation of which is standard. We use the signature and sign conventions of 0. The 
sum in ([!]) is taken over the bulk components bounded by branes, and Aj is the cosmological 
constant in the zth bulk component. Note that we have included the curvature term in 



the action for the brane with the coefficient m 2 . Such a term generically arises when one 
considers quantum effects generated by matter fields residing on the brane in the 

spirit of the idea of the induced effective action [H||; its inclusion in the action is therefore 
mandatory. The lagrangian L (h a p, <p) corresponds to the presence of matter fields <f> on the 
brane interacting with the induced metric h a p and describes their dynamics. The extrinsic 
curvature K a g on both sides of the brane is defined with respect to the inner normal n a , as 



m 



Let us first consider the case with Z2 symmetry of reflection, which requires equal cos- 
mological constants on the two sides of the brane: Ai = A2 = A. The corresponding 
cosmological equation of theory ([!]) was derived in and has the form 



a 2 3m 2 I 2 



1± \ 1 



p + cr A C 



(2) 



3m 2 6 a 4 , 

where the integration constant C corresponds to the presence of a black hole in the five- 
dimensional bulk solution, and the term C/a 4 (occasionally referred to as 'dark radiation') 
arises due to the projection of the bulk gravitational degrees of freedom onto the brane. The 
length scale t is defined as 

< 3 » 

The new singularities that we are going to discuss in this paper are connected with the 
fact that the expression under the square root of (@) turns to zero at some point during 
evolution, so that solutions of the cosmological equations cannot be continued beyond this 
point. There are essentially two types of 'quiescent' singularities displaying this behaviour: 

A type 1 singularity (SI) is essentially induced by the presence of the 'dark radiation' 
term under the square root of @ and arises in either of the following two cases: 

• C > and the density of matter increases slower than a -4 as a — > 0. Such singularities 
occur if the universe is filled with matter having equation of state P/p < 1/3, an 
example is provided by pressureless matter (dust) for which p cx a~ 3 . A special case 
is an empty universe (p = 0). 



• The energy density of the universe is radiation-dominated so that p = po/a 4 and 
C > p . 



The singularities discussed above can take place either in the past of an expanding 
universe or in the future of a collapsing one. 
A type 2 singularity (S2) arises if 



In this case, it is important to note that the combination p/3m 2 — C/a A decreases mono- 
tonically as the universe expands. The expression under the square root of (0) can therefore 
become zero at suitably late times, in which case the cosmological solution cannot be ex- 
tended beyond this time. S2 is even more interesting than SI since: (i) it can occur during 
the late time expansion of the universe; (ii) it can occur even if dark radiation is entirely 
absent (C = 0). 

For both SI & S2, the scale factor a(t) and its first time derivative remain finite, while 
all the higher time derivatives of a(t) tend to infinity as the singularity is approached. As 
an example consider a type 2 singularity with C — 0, for which 



as p(t) — > a = Am 2 /2 — o — 3m 2 ji 2 . We therefore find that the scalar curvature R — > oo 
near the singularity, while the energy density and pressure remain finite. Although this 
situation is quite unusual from the viewpoint of the intrinsic dynamics on the brane, it 
becomes comprehensible when one considers the embedding of the brane in the bulk. It 
is well known that the cosmological braneworld under consideration can be isometrically 
embedded in the five-dimensional solution of the vacuum Einstein equations described by 
the metric 




(4) 




n > 2, 



(5) 




(6) 
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where 

/»(0 = (7) 

and dVL K is the metric of the three-dimensional Euclidean space, sphere or pseudosphere 
corresponding to the value of k = 0, ±1. The embedding of the brane is defined by the 
function r = a(t), and one can then proceed to define evolution in terms of the proper 
cosmological time r of the induced metric on the brane. The cosmological singularity under 
consideration is connected with the fact that the brane embedding is not extendable beyond 
some moment of time t because the function a(t) that defines the embedding cannot be 
smoothly continued beyond this point (see Fig. |I|). 




FIG. 1. Conformal diagram showing the trajectory of a spatially spherical braneworld embed- 
ded in the five-dimensional Schwarzschild space-time. The trajectory is not smoothly extendable 
beyond the point Q. 




FIG. 2. The involute 2 of a planar curve C is not smoothly extendable beyond the starting 
point Q. 
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This specific feature of the brane embedding can be compared to the behaviour of the 
involute of a planar curve. The involute X of a convex planar curve C is a line which intersects 
the tangent lines of C orthogonally ||. X can be visualised as the trajectory described by 
the end of a strained thread winding up from C (see Fig. |^). The involute of a typical curve 
is sharp at the starting point Q so that it is not smoothly extendable beyond the point Q. 

This analogy can be traced further. Note that the evolution of the brane in theory ([I]) 
is described by the following well-known equation: 

m 2 G af3 + ah a p = T aP + M 3 ^ {K aP - h aP K) , (8) 

where the extrinsic curvature is summed over both sides of the brane. One can see that 
it is the influence of the sum of the extrinsic curvatures on the right-hand side that leads 
to the singularities under investigation so that the singularity of the Einstein tensor G a p is 
accompanied by the singularity of the extrinsic curvature K a p, while the induced metric h a p 
and the stress-energy tensor T a p on the brane remain finite. Quite similarly, the involute of 
a curve is defined through the extrinsic curvature of its embedding in the plane, as is clear 
from Fig. 0, and its singularity at the point Q is connected with the fact that the extrinsic 
curvature diverges at this point. Specifically, the parametric equation for the involute x*(s), 
s > 0, in Cartesian coordinates on the plane can be written as [[| 

x*(s) = x(g — s) + s ■ x'(g - s) , (9) 

where x(s) is the curve C parametrized by the natural parameter s, and x(g) = x*(0) is the 
coordinate of the starting point Q of the involute. The extrinsic curvature of the involute is 

k(s) = -, (10) 
s 

which diverges at the starting point Q corresponding to s = 0. 

One should also highlight an important difference between the ID and 4D embeddings: 
the involute being one-dimensional, a singularity in its extrinsic curvature does not lead to 
a singularity in its intrinsic geometry. As we have seen, this is not the case with the brane 
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for which the extrinsic and intrinsic curvatures are related through ([§]) , so that a singularity 
in K a/3 is reflected in a singularity in G a @. 

Interestingly, an S2 singularity can arise in the distant future of a universe resembling 
our own ! To illustrate this we rewrite Eq. (0) (with C = 0) as || 

H\z) 



H 2 



n m (i+ z y + n K (i+ z y + + 2n e ± 2Jn e Jn m (i+ Z ) 3 + n a + n e + n A , (11) 



The underlining emphasises those terms which cause the braneworld to differ from its 
general- relativistic counterpart. For simplicity, we shall only discuss the solution corre- 
sponding to the '+' sign in (|TTD (called BRANE2 in ||). In this case, it can be shown that 
our model satisfies the constraint equation 



n m + n K + n a + 2yjn e yji-n K + n A = i, (12) 

where 

O ... Po n «_ n ° n - 1 n ^_ (io\ 

m ~3m 2 tf 2 ' alHf CT ~ 3m 2 tf 2 ' l ~PH*' 6# 2 ' 1 J 

and the subscript '0' refers to the present value of the various cosmological quantities. 

Inequality (Q) now becomes 

n a + n e + n A < o , (14) 

and the limiting redshift, z s = ao/a(z s ) — 1, at which the braneworld becomes singular is 
given by 



The time of occurance of the singularity (measured from the present moment) can easily be 



determined from 



r° dz 

r - = «<* = *>-'<* = °> = 1 (16) 



where H(z) is given by ( |TTj ) (see also [IU|). In Fig. [5] we show a specific braneworld model 
having Q m = 0.2, = 0.4, Q A = Q K = 0. In keeping with observations of high redshift 
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0_ = 0.2 



q(z) 



redshift z 

FIG. 3. The deceleration parameter (solid line) is shown for a braneworld model with S7 m = 0.2, 
= 0.4, = = 0, and VL a determined from (|l^). We find that q(z) — > 0.5 for z 3> 1 while 

5(2) — > 00 as z — > —0.312779... Currently qo < 0, which indicates that the universe is accelerating. 

Also shown is the dimensionless Hubble parameter h(z) = 0.1 x H(z)/Hq (dashed line) for this 

model. The vertical line at z = shows the present epoch. 



supernovae our model universe is currently accelerating ||11|| , but will become singular at 
z s ~ —0.3 =>- a(z s ) ~ 1.4 x a , i.e. after T s ~ 4.5 Gyr (h wo = H o /100 km/sec/Mpc). 
Figure demonstrates that the deceleration parameter becomes singular as z s is approached: 



d H' 



(17) 



while the Hubble parameter remains finite: 

H 2 (z.) 



H 2 



A • 



(18) 



It should be noted that, for a subset of parameter values, inequality (|14]) can be satisfied 
simultaneously with > ■ In these models, the universe will recollapse (under the 
influence of the negative brane tension) before the S2 singularity is reached. (The marginal 
case Q\ = corresponds to the Hubble parameter vanishing at the singularity.) 
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Now let us briefly consider the general braneworld without Z<i symmetry. In this case, 
theory (TT]) leads to the following general cosmological equation for the brane embedded into 
the five-dimensional bulk [HI: 



M Ul K P+*V _ AA/ r6 (tt2 , K Ax+A 2 C 



m ( H + — r — V-s- = 4M b iZ* + 



a 2 3m 2 / V a 2 12 



M 



12 



36m 4 



Ax - A 2 + E/a 



i 



n 2 



H 2 + K/a 2 -(p + a)/3m 2 



(19) 



where C has the same meaning as in (g), is another arbitrary integration constant, and 
Ai and A 2 are the cosmological constants on the two sides of the brane. In this paper, we 
restrict ourselves to the situation where Ai = A 2 = A, but E ^ 0. One should note that 
(H) reduces to <g) if A x = A 2 = A and E = g§]. 

It is easy to see that singularity S2 is always present in the past of the expanding brane. 
The reason for this rests in the negative character of the last term on the right-hand side of 
Eq. (|T9|), which rapidly grows by absolute value as a —> 0, while the left-hand side of this 
equation is constrained to remain positive. In the case of an expanding brane (assuming 
Ai = A 2 = A), the last term on the right-hand side of (|T9|) rapidly decays and becomes 
unimportant. Therefore, provided (f|) is satisfied, the expanding universe will encounter an 
S2 singularity in its future. 

We note that all the singularities described above are absent in the limit m — > which 
is frequently discussed in the literature. Indeed, in this limit, Eq. (|l!|) takes the form 



H + ^ = ^2— + ^ + ^6M^ + T^l -pTa 1 ' (20) 



which only admits cosmological singularities associated with an infinite density of matter 
and dark radiation (C/a 4 ). These singularities are reached as a — > 0. (In general relativity, 
homogeneous and isotropic space-times generically admit only infinite-density singularities, 



whereas anisotropic space-times can be empty yet singular ||12|| .) 

Singularites of the kind discussed in this paper occur when the original equations of 
motion are non-linear with respect to the highest derivative. They have earlier been discussed 



in the context of Einstein gravity with the conformal anomaly ||13|| . (This result is not 
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surprising in view of the formal similarity between braneworld theory and GR-based models 
with the conformal anomaly, discussed in ||.) 'Determinant singularities' having a similar 
structure and properties are known to arise in the anisotropic Bianchi I model containing 
a dilaton coupled to a Gauss-Bonnet term in the action [[TJj] (the general theory of such 
'determinant singularities' is discussed in JT3|). We should emphasise that both the past 
(SI) and future (SI, S2) singularities in the braneworld scenario occur for a wide range in 
parameter space and might provide an interesting alternative to the 'big bang'/'big crunch' 
singularities of general relativity. In addition, the universe of Fig. |] which terminates at 
the S2 singularity with p m — > constant as t — > t s , provides an interesting contrast to the 
bleak finale presented by quintessence or cosmological-constant dominated models, in which 
Pm —> as t — > oo. 

Finally, we would like to draw the reader's attention to issues which we feel require further 
investigation. First of all, the singularities discussed in this paper occur for the very special 
class of homogeneous and isotropic spaces. The singularities in this case are characterized by 
the property that the energy density remains bounded while the curvature blows up as one 
approaches the space-time singularity. It is open to question whether singularities of this kind 
will persist in space-times having less symmetry, including the more general anisotropic and 
inhomogeneous case. Secondly, as the singularity is approached and the scalar curvature of 
the induced metric on the brane diverges, higher-order derivative terms in the effective action 
for the brane become important. Taking these into account may qualitatively modify the 
behaviour of solutions as is the case in general relativity. Thirdly, back-reaction associated 
with quantum effects such as particle production and vacuum polarization may influence 
the behaviour of cosmological solutions in the vicinty of the singularity. We hope to return 
to these issues in our subsequent work. 
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